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Comment on "Disorder Induced Quantum 
Phase Transition in Random-Exchange 
Spin-1/2 Chains" 

In a recent Letter Hamacher, Stolze and Wenzel [1] 
studied the disordered spin-1/2 AFM XXZ chain 


H = J 


L-l 


^i(SfSf +1 + SfSf +1 ) + AS*Sf +1 


with Xi i.i.d. random variables uniformly distributed over 
the interval [1 — W; 1 + W], the parameter W controlling 
the strength of the disorder. The authors claim that for 
A < 1 they found numerical evidence for non-universal 
behavior for weak disorder (W < 1) manifested in a con- 
tinously varying exponent r)(W) describing the asymp- 
totic decay of the transverse spin correlations 


around L/£w ~ T the other estimates for £w are then 
chosen to give the best data collapse. We see that for all 
disorder strength the maximum system sizes used in [1] 
are still well within the crossover region and not in the 
asymptotic regime. We do not expect that this situation 
will improve for A > 0. In general £w diverges when 
W — ► 0, and we found that our estimates for £vv obey 
(1) Zw 

oc where S is the (VF-dependent) variance of the 
random variable In A; , which determines according to [2] 
the RG-flow of the model (1). We obtain $ « 1.24±0.02. 


C xx (r) = (Sf 5f +r ) oc r~" (M/) 


(2) 


where (...) denotes the ground state expectation value 
averaged over the disorder and the sites i. They con- 
cluded that there is no universal infinite randomness fixed 
point (IRFP) as predicted by D. Fisher [2]. 

In this comment we show that these conclusions are 
wrong and that the numerical data presented in [1] are 
not in the asymptotic regime. This misconception is due 
to the fact that the authors ignored the existence of a In- 
dependent crossover length scale £w that describes the 
crossover from the pure fixed point to the only relevant 
IRFP: For L <C £w one observes the critical behavior of 
the pure system (Aj = const.), and only for L 3> £,w the 
true asymptotic critical behavior (r)(W) = 2, indepen- 
dent of W > 0) of the disordered chain becomes visible. 
Even for strong disorder (W = 1.0) £w is of the same 
order of magnitude as the system sizes considered in [1]. 

In order to be able to reach sufficiently large system 
sizes we restrict ourselves to A = in which case (1) 
reduces to a free fermion model and the ground state 
computations are done following [3]. In Fig. 1 we show 
the averaged bulk correlation function C xx (L/2) for dif- 
ferent strengths of the disorder. We observe that asymp- 
totically (i.e. for L — > oo) the data follow the behavior 
C XX (L) oc L~ 2 as predicted by the real space renormal- 
ization group [2]. Only for small L the data appear to 
follow a non-universal (i.e. VF-dependent) power law, and 
this is the region on which [1] reports. 

Due to the presence of a crossover length scale £(W) 
the correlation function obeys the scaling form 


C°*(L/2) = L-^ciL/tw) 


(3) 


where c(x)const for x — > 0, and c(x) — > x 3 / 2 for x — ► oo. 
This implies C xx {L/2) oc L^ 1 / 2 for L < £ w (the pure 
behavior) and C xx {L/2) oc L~ 2 for L > £ w (the IRFP 
behavior). In the inset of Fig. 1 we show such a scal- 
ing plot of the data in the main figure. We have chosen 
Zw for W = 1 such that the crossover region is centered 


0.1 


0.01 


0.001 


0.0001 


1e-05 


W = 1 — 

W = 0.25 

W = 0.50 □ 

W = 0.75 — ■- 

w = 1 .q - - e - 


0.01 



c w - 600 
SJ-140 


10 


100 
L 


1000 


10000 


FIG. 1. Averaged correlation function C xx (L/2) as a func- 
tion of the system size L in a log-log-plot for W=0, 0.25, 0.5, 
0.75, 1.0 (top to bottom). The data are averaged over 50000 
for L < 1024, 3000 for L=2048 and 500 for L=4096, the statis- 
tical error is smaller than the symbol sizes. The data for the 
pure system {W = 0) follow C xx (L/2) oc L" 1/2 , the full line 
with slope —2 is the expected asymptotic behavior according 
to the IRFP scenario [2]. Inset: Scaling plot according to 
Eq. (3) for the data of the main figure with £w ~ 600, 140, 
54, 20 for W = 0.25, 0.5, 0.75 and 1.0, respectively 

To conclude we have shown that contrary to what has 
claimed in [1] there is no non-universal behavior in model 
(1) for weak disorder and that by taking into account the 
existence of a crossover length scale the numerical data, 
ours as well as theirs, for A < 1 are compatible with the 
IRFP scenario predicted in [2]. 
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